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Abstract 

Deformations that conserve the parallehsm and the distances -between layers, 
in smectic phases; between columns, in columnar phases- are commonplace in 
liquid crystals. The resulting deformed textures have the same mass density 
as in the ground state (an expected property in a liquid) and are at the same 
time isogonic and isometric, which imposes specific geometric features. The 
corresponding order parameter singularities extend over rather large, macro- 
scopic, distances, e.g., cofocal conies in smectics. This well-known picture is 
modified when, superimposed to the ID or 2D periodicities, the structure is 
helical. Isogony is no longer the rule, but isometry (and mass density) can be 
preserved. This paper discusses the case of a medium whose structure is made 
of ID modulated layers (a lamello-columnar phase), assuming that the mod- 
ulations rotate helically from one layer to the next. The price to pay is that 
any isometric texture is necessarily frustrated; it consists of layers folded into 
a set of parallel helicoids, in the manner of a screw dislocation (of macroscopic 
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Burgers vector), the modulations being along the helices, i.e. double-twisted. 
The singularity set is made of two helical disclination lines. We complete this 
geometric analysis by a crude calculation of the energy of a helical ribbon. It 
is suggested that the helical ribbons observed in the B7 phase of banana-like 
molecules are such isometric textures. As a side result, let us mention that 
the description of double- twist, traditionally made in terms of a partition of 
the director field into nested cylinders, could more than often be profitably 
tested against a partition into nested helicoids. 
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So full of shapes is fancy, 

That it alone is high fantastical. 

Twelfth Night [I, 1, lines 14-15] 

I. INTRODUCTION 

Consider a bent solid crystal (Fig. la); it is well known since [1] that the bending 
stresses can be relaxed by a set of dislocations of density 1/bR, where R is the radius of 
curvature and b the Burgers vector. These dislocations (in this example, edge dislocations) 
are spaced at a distance ~ VbR. Under annealing they assemble in the shape of walls; 
this process is known by the name of polygonization [2,3]. In this geometry (Fig. lb) the 
layers of atoms parallel to the Burgers vector keep practically unstrained, and the largest 
contribution to the strain energy comes from the cores of the dislocations. The repeat 
distance in a direction perpendicular to the Burgers vector is conserved and the atoms have 
almost the same density along the curved layer as in the unbent crystal; the curvature can 
be large. The same geometry transposed to a layered liquid, a smectic A (SmA) , shows that 
the curvature deformation is compensated by liquid relaxation along the layers. Fig. Ic [4] . 
Transposed to a columnar phase, the columns curvature is such that their equidistance and 
parallelism are conserved, which implies that they are orthogonal to a family of planes; the 
2D network of their intersections with any of those planes is unchanged compared to the 
undeformed state [4]. In a sense it is still feasible to speak of dislocations in lamellar and 
columnar phases, but they now have a vanishing Burgers vector, their density is infinite, and 
there is no analogue to polygonization. The energy of deformation is pure curvature energy; 
the mass density is conserved. Such deformations, we shall call them isometric, in the sense 
that they do not affect the 'sohd' component of the ordered hquid, i.e. the distance between 
lamellae or between columns. In the columnar case, they also are isogonic. 

Fig.l 
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Isometric deformations are commonplace in ID and 2D liquid crystals. It is an easy 
matter to show that they must dominate on scales larger than \J^-, where X is a typical 
Frank curvature modulus and B a typical inverse compressibility, i.e. on scales larger than 
a few repeat distances [4]. They form characteristic textures that carry singularity sets. 
The case of smectics has been known for more than a century [5]; the relevant isometric 
textures are focal conic domains., which are singular on a pair of conjugate conies, the layers 
take the shape of Dupin Cyclides. In columnar phases the isometric textures are developable 
domains [6,7], which are singular on a developable surface; the columns can be partitioned 
(in an infinity of ways) into families of parallel lines belonging to Monge Surfaces. We refer 
the reader to Ref. [4], chapter 10, for a full description. It is a remarkable fact that these 
objects, whose geometry is so precisely defined, extend over mesoscopic scales, quite often of 
the order of a mm, i.e. involving 10^ to 10^ strictly parallel layers in a SmA. This indicates 
the extension of the isometric principle in ordered liquids. 

Let us come now to a newcomer in the realm of mesomorphic structures, namely the 
B7 phase. It is widely accepted that the structure of this phase involves at the same time 
a smectic layering (Sm) and a one-dimensional modulation (D) of the layers [8] (we call 
'columns' the elements of this modulation, in a metaphoric sense; we denote by it their 
unit tangent vector; the symbol D stands for 'order of the columnar type', as with discotic 
thermotropic molecules). It is believed that the B7 order is built upon a low-symmetry, 
chiral, molecular configuration [9], although the molecules are non-chiral. These molecules 
arc planar and bent-like, with a rather strong dipole moment along an axis of symmetry. 
Fig. 2. Therefore domains of opposite chiralities coexist in the ground state. But at a 
mesoscopic level, notwithstanding this chiral character of the structure, the ground state 
does not exhibit helicity (only Sm and D orders are simultaneously present). Observe that 
the lamellar structure allows that the modulations of neighboring layers keep in register, 
yielding thus close-packing and high mass density. 

Were it limited to these two types of order, Sm+D, one would expect that the textures 
of lower energy, i.e., the isometric distortions, although somewhat more subtle than the 

4 



isometric distortions of sets of parallel 2D liquid layers or of parallel ID liquid columns, 
would not require an analysis of a degree of difficulty much deeper, in the sense that no new 
geometrical concept would be necessary [10]. 

A most conspicuous feature of the B7 phase is the presence of helical ribbons (HR) that 
appear in the isotropic phase when the samples are cooled down rapidly from the high 
temperature range of the isotropic phase [8,11-14]. These nuclei clearly display a double- 
helix texture. Fig. 3. Left and right HRs are in approximately equal numbers. A linear 
relation between the pitch b and the width I has been documented by several authors [8,13]. 
The model for HRs proposed by Coleman & al. [14] consists in two cylinders, made of tightly 
packed cylindrical layers, twisted one about the other. This model is in full agreement with 
the Sm+D structure of the B7 phase, but is not isometric, in the sense that the twisting of 
the two cylinders hinders close-packing. 

Fig, 2 

There are in fact two types of HRs [12]. The S-HRs (slim) are characterized by a rather 
perfect double hehcal geometry, a pitch/width b/l ratio ^ n (see Fig. 3), and a high rigidity; 
they do not bend easily, which is related to the remarkable constancy of the pitch and of 
the b/l ratio all along each ribbon. The F-HRs (fat) are much more disperse in pitch, shape, 
and width. We shall in this note expatiate only on S-HRs, (simply noted HRs in the sequel). 
We present in the next section a purely geometric model, which interprets HRs as isometric 
textures of the B7 phase; and explains satisfactorily well their double helical geometry and 
the value of the ratio 6/Z = tt. In particular, the model is such that the modulation of any 
B7 twisted layer does match the modulation of the neighboring twisted, parallel, B7 layers, 
as required by the isometric principle. 

Fig.3 

The helical pitch b varies somewhat from one nucleus to another, but about a value 
that is rather strongly peaked. Fig. 4. Also the S-HRs transform in a period of a few 
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minutes into F-HRs. These features have led the authors of [12] to assume that S-HRs 
are specific textures of metastable modifications of the B7 phase, which display helicity (H 
order), either left or right, superimposed on Sm and D orders. We note them B7*. The 
axis of helicity is perpendicular to the layers and affects the columns. We thus assume that 
the B7* hehcal states differ from B7 by the existence of a finite pitch p of the columns, p 
being a material constant. We show, in the framework of a very simple model, that the free 
energy is minimized when voq = —1 (where vo — 6/27r, q — 27r/p), i.e. when h = — p, which 
is precisely what is expected from the purely geometric model, and that this energy is less 
than the energy of the model advocated in [14]. 

Fig.4 

The next section is devoted to a detailed account of our HR geometric model; section 
III deals with an approached calculation of the HR energy, assuming the existence of a 
metastable B7* phase. 

II. THE GEOMETRICAL MODEL. 

A. Giant Screw Dislocations in SmAs, the Chromosome of Dinoflagellate, 

Double- Twist, and all that. 

Our isometric model is related in many ways to the three topics mentioned in the title 
of this subsection. We briefly expose this relationship with these topics, before (in the next 
subsection) giving a more detailed account of it. 

Let us first summarize the main features which characterize a HR according to our model: 
a HR is the inner part of a giant Burgers vector screw dislocation; the layers are deployed 
along parallel and equidistant helicoids whose common axis is the dislocation line and whose 
pitch is equal to the Burgers vector; the columns are along parallel hehces of the same axis 
and belong to the helicoids. The singularity set of this isometric texture is made of a double 



6 



helix located at the periphery of the whole texture. Fig. 5 is a mathematical illustration of 
the model; see also Fig. 8 of [12]. 



Screw dislocations with giant Burgers vectors are documented in SmAs since long [15]; 
they are characterized by the fact that the core is split into two helical singularities located 
at a distance w — hjli: from the axis of the screw, h being its Burgers vector; this geometry 
is visible under the optical microscope, h being so large (typically a few 10//m in CBOOA, a 
classic thermotropic smectic; the Burgers vector is not related to any material constant, but 
rather to the circumstances of the experiment, very probably the thickness of the sample). 
This geometry can be explained as follows [16,4]. The SmA layers form a set of parallel, 
equidistant smf aces (isometric principle!), which are nested helicoids, stacked upon a central 
ruled helicoid^ whose axis is along the screw line; it has a right or left hand, according to 
the case, analytically distinguished by the sign of b. A ruled helicoid is generated by a 
straight-line meeting a fixed axis at a right angle and moving with a uniform helical motion 
about this axis; it has a vanishing mean curvature: 



The two singularities are the cuspidal edges of the focal surfaces of the set of helicoids; 
they act physically as two helical disclinations of strength k — 1/2 (see [4]), about which 
the layers fold. Outside the cylindrical region bound by the disclinations, the layers adopt a 
quasi-planar geometry that is not in analytical continuity with the helicoids inside, but are 
in register with the layers inside. 



^It does not come as a surprise that the same object is sometimes called 'right helicoid' in Anglo- 
Saxon countries (see for example [17]), and 'helicoi'de gauche a plan directeur' by the French [18]; 
of course, in these expressions, neither the attributive 'right', nor the epithete ^gauche' refer to 
handedness! We use a third denomination. 



Fig. 5 




(1) 
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We argue that the B7* layers in the HR stack as those of the SmA phase in the screw- 
dislocation. The geometric properties of the layer stacking impose that the disclination lines 
are scribed on a cylinder of diameter I = b/ir. This is precisely what one observes empirically 
in the B7* case, taking for I the measured diameter of the HR. But there are differences. The 
HR is limited to the core region, inside the cylinder on which the disclinations are scribed, 
whereas the SmA screw dislocation extends outside the core region. Also, it is precisely the 
hehcity of the columns (no columns in the SmA phase, of course) that stabilizes the HR. 
The screw dislocation is a metastable defect of the SmA phase. The HR, on the other hand, 
see next section, is the texture of lowest energy of the frustrated metastable B7* phase. 

The HRs bear also some resemblance with a model [19] devised in order to explain the 
conformation of the molecule of DNA in the chromosome of Dinofiagellate, observed by 
Livolant & Bouligand [20]. The same observations also attracted the attention of J. Friedel, 
who concluded to a much similar model [21]. As shown in [20], the DNA is cholesteric in this 
chromosome, but at the same time the isometric principle, which applies here because the 
persistence length of the molecule is not negligible compared to the cholesteric pitch, requires 
that the local packing be as dense as possible. Ref [19] discusses precisely the possibility 
of the coexistence of cholesteric and 2D local order, which a priori are incompatible in the 
large, but, as it was shown, can be fitted locally if the molecules can be partitioned into 
a family of parallel, helical, surfaces containing their long direction. Here sits the analogy 
with HR. But whereas the columns are along hehces in a HR, it is beheved that the DNA 
molecules (but to shape fluctuations) are along a set of lines orthogonal to these helices, 
namely a set of geodesies^. 

In B7* the columns sit along helices with approximately the opposite pitch to the layers 
pitch; thus, they display the well-known blue phase (BP) double-twist geometry. In the case 



^This is not the densest packing; but it is the straightest configuration, which is in accord with 
the long persistence length of DNA. 
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of BPs, this geometry is not related to the presence of any layers: the columns are in fact 
liquid-like in the three directions, and the order parameter is locally cholesteric. Double- 
twist is splayless and its bend energy is relatively small, as long as the lateral size of the 
bundle of double-twisted molecules does not exceed ^ b/2 (the double- twist geometry is 
frustrated). It has been stressed [22,23] that in ideal double- twist geometry, the columns are 
equidistant, having the same property of equidistance met e.g., for the strands of a string. 
It is such a property we expect that the columns of the B7* phase do satisfy, with the result 
that the elastic energy is not prohibitive. How this requirement relates to helical layers is 
discussed now. But whereas most people think of the double-twist geometry as geometry of 
nested cylinders, into which the columns are partitioned, the screw dislocation model drives 
us to partition the HR geometry into hehcal surfaces. We shall come again to the nature of 
this discrepancy, which is not as trivial as it might look at first sight. 

B. Why Helicoids 

Helicoids possess a property that they share with surfaces of revolution: they are invariant 
under a one-parameter family of rigid motions [17], namely a continuous set of translation- 
rotations about the helical axis (for helicoids), a continuous set of rotations about the axis 
of revolution (for surfaces of revolution). Of course the circular cylinder (limit case of a 
helicoid for w = 0) and the sphere possess respectively two- and three-parameter families of 
rigid motions, but we leave aside these limit cases, that are trivial. Less trivial is the fact 
that surfaces of the type under consideration, namely helicoids (resp. surfaces of revolution), 
contain two orthogonal families of curves of which one family is the family of helices (resp. 
circles of latitude) whereas the other one consists of geodesies (resp. meridians). This is 
exactly the geometrical property we are looking for: if one installs the columns of the B7* 
phase along the curves orthogonal to the geodesies, then the columns stay at equal distances, 
because of the property of frontality of the geodesies. 

The demonstration of the existence of two orthogonal families of curves [17] requires 
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the recognition that the rigid motion invariance is also invariance under a one-parameter 
family of bendings. A surface admits bending invariance if the following is true. Place, snug 
against the surface, a piece of flexible but not stretchable foil; then it can be displaced along 
the surface, at least along one direction, so that it keeps snug against it without tearing, 
although it can change its shape. If the surface admits a one-parameter family of bendings, 
then there are curves, which cover the whole surface, such that the displacements along those 
curves exhibit the above property of invariance. Because bending does not change geodesic 
distance, any two curves of the family must be at a constant geodesic distance. Therefore a 
geodesic line orthogonal to one of the curve of the first family crosses at right angles all the 
curves of the first family. Because bending keeps Gaussian curvature invariant, Gaussian 
curvature is constant all along a curve of the first family. Gaussian curvatures measured 
along two curves of the first family are generically different. Therefore the curves of the first 
family do not cross; they form a family without singular points. 

This first family of lines is a set of helices. It can indeed be shown [17] that surfaces 
admitting a one-parameter family of bendings can be brought into the form of general 
hehcoids or of surfaces of revolution. Most evidently the curves of the first family are 
the helices (resp. the circles of latitude), and no other curves. Therefore the condition of 
equidistance (more precisely: of equal geodesic distance) forces the columns to sit along 
the hehces of general helicoids (resp. along the circles of latitude of surfaces of revolution). 
The physical layers, which are deployed along these helicoids, have to be parallel in order 
to minimize compression energy. We are thus led to study sets of parallel general helicoids, 
and their helices. 

The property of columns equidistance, which is preserved in each individual helicoid, 
is still true for columns belonging to different parallel helicoids. In effect, consider two 
neighboring parallel surfaces H and H' at a distance SX measured along their common 
normals. Let 

(71(72, ^((7l + (72) (2) 
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be the Gaussian and mean curvatures at some point M on H. G and H are constant along 
the hehx h passing through M on H, and depend upon the unique parameter that defines the 
hehx. Two parallel surfaces have the same centers of curvature Ci and C2 at corresponding 
points, say M (on H) and M' (on H'). Therefore 

" 1 + 5Xa, ' '^^ " 1 + 5Xa2 ^ ^ 

Thus the Gaussian and mean curvatures G' and i?' are constant along the curve h' drawn 
on the surface H' by the normals leaning on the first helix h. Therefore H' is invariant under 
a one-parameter family of bendings; it is indeed a helicoid, as already stated. On the other 
hand, h', which by construction is at a constant distance SX of h, is an helix on H', and thus 
belongs to a family of equidistant curves on H', all obtained by tracing the curves h'(5A) 
equidistant to the helices h. 

This property explains why the isometric principle is satisfied not only on each helicoid 
separately, but also throughout the whole HR. Notice that the columns are not parallel, but 
anyway equidistant, which allows the modulations to fit tightly from one layer to the next. 
It is fully evident that a tight fitting is not feasible with planar parallel layers. 



C. Analytical Description of a Simple Stacking of Helicoids, Ingredients 

A parallel stacking of helicoids on a ruled helicoid has been considered in [19]. We revisit 
the corresponding analytical description, but the detailed calculations, already in [19], will 
not be given. The next section will be devoted to the calculation of the energy of such an 
assembly of layers, and its comparison with a cylindrical stacking, which yields a somewhat 
comparable but in fact fundamentally different assembly of helical columns. 

This subsection gives the analytical expressions necessary to calculate the energy of a HR 
(next section); it is mostly technical, and can be skipped in a first reading. Note however 
the equations 4 and 5, for the notations which are used; they give the expressions of the 
central ruled hehcoid H(0) and of any parallel helicoid H(A) at a distance A. Equations 9 



11 



and 10 give the coordinates of the tangent unit vector it along the hehces (the columns); 
the important point is that it is rotationally invariant about the z-axis, which lets suppose 
that the description of a HR in terms of nested cylinders is as valid a description as in terms 
of nested helicoids. In fact, it is not so; we will discuss later on why. 

Let A be the signed distance of the hehcoid H(A) to the central right helicoid H(0). A 
point M on H(A) is identified by the r and 9 values of the foot on H(0) of the normal to 
H(A) in M. Because the hehcoids are parallel, they all have the same set of normals. The 
coordinates of a point on H(0) are 

H{0) = {x ^ rcos9,y ^ rsm9,z ^ 'nj9}; (4) 

those of a point on H(A) are 

H{X) = {x ~ rcos9 + ^^sm9,y — rsm9 — ^^cos9,z — w9 + ^}, (5) 
N"^ ^r'^ + w'^ . (6) 

We introduce the notation t = N'^/Xw. The Gaussian curvature and the mean curvature 
at M can be written 

_ _ _ 

^ - - ~it^-i)m - ~x^{t^-iy 

One observes that: 

(1) - the mean curvature vanishes for A = 0, i.e. on the central helicoid H(0), Fig. 5b. 

H(0) is indeed a (ruled) minimal surface. Note that the sign of the mean curvature depends 

on the orientation of the normal to H(A); here we have assumed that the three vectors -r— , 

or 

-— -, and the normal u to H(A) in M form a right-handed frame of reference. 
d9 

(2) -the Gaussian curvature is negative for {t^ — 1) > 0, positive otherwise, and is infinite 
for t = ±1. The region of negative Gaussian curvature contains the central hehcoid - since 
it has negative Gaussian curvature {t = ±oo) - and extends symmetrically on both sides of 



12 



H(0) until it meets the two focal surfaces {t — ±1) of the family of hehcoids, i.e., the envelop 
of their normals. These focal surfaces were studied in some details in [19]; see also [4]; we 
recall that each of them possesses a helical singular line, around which the layers show a 
k — +1/2 configuration. 

(3)-G, H, and ai, a-i depend only on r (not on ff) i.e. on the parameter that is a constant 
on a helix. Therefore the curvatures properties, in particular the Gaussian curvature, are 
constant along a given helix, which agrees with our analysis above. 

The coordinates of the tangent unit vector along the helices can be written: 

it = {cos 6 cos V ~ ^ ^ cos V cos 9 sin (/?, — sin (/?}; (9) 

(/? given hj t — tanc/? = N^/Xzu. 

Using Eq.5, ft can also be written 

It = ^^{-y,x,zu}. (10) 

Hence it is a splayless {divlt = 0) double-twisted director field, rotationally symmetric 
about the 2;-axis. The expressions below will be used in evaluating the energy of the system: 

it ■ curllt = —2(p', \lt X curlltf — (p'^ — (p'^, (11) 

it X curllt — — simp {uy, —Ux, 0} , \curlltf — (p'^ + 3</?^, 

where ^'\ — sinip). The first equation on the second line requires — 1 < vop)')^ < 0. In 
the applications, we shall assume, without loss of generality, cc/ > 0, ip'^ < 0. Thus sin 2ip > 0, 
and < if < 7r/2, but because |tan(^| > 1, one gets: 7r/4 < cp < 7r/2. Furthermore, because 
of the symmetry A —A, the points M+ on H(-|-A) and M_ on H(-A) are coupled by pairs 
that are symmetric with respect to the intersection of the segment M_|_M_with H(0) in its 
middle; we restrict in the sequel to the region A > 0. 

There are two types of hehcoids H(A), those that are everywhere of negative Gaussian 
curvature {X < zu), those whose Gaussian curvature changes sign [X > zu). The sign of the 
Gaussian curvature is negative ii > 1, which can also be written > —zu'^ + Xzu. 
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(1) - Helicoids X < zu 

< r < oo; arctan ^ <(p < ^. (12) 

These helicoids do not touch the focal surfaces. They are folded around the 
hehcal dischnation k = 1/2 that corresponds to the helix ip — 7r/4, r = 0, 
{x = w sin 9,y = —zu cos 9,z = —w6} on the helicoid H(A — w). The Gaussian curvature is 
negative over the entire helicoid, when \ < vo. 

(2) - Helicoids X> zu 

One distinguishes two regions 

(a) - negative Gaussian curvature: 

zu {X — zu) < r"^ < oo; — <(fi < —. (13) 

(b) - positive Gaussian curvature 

Q <r^ <w{X — w)] arctan ^ <ip < j. (14) 



D. Why Negative Gaussian Curvature 

The helicoids are never physically present in their entirety in a HR; and the physical 
parts of the helicoids are those which have negative Gaussian curvature. The regions with 
positive Gaussian curvature are disfavored, and the span of the helicoids is finite, for the 
following reasons: 

(a)- The helicoids X < zu (which have G < 0) are free of any singularity. The helical 
dischnation k — 1/2 , which is the first singularity met when proceeding away from the 
central axis, belongs to the helicoids X — zu. But the helicoids X < zu do not extend very 
far from the central axis because, as in double-twisted cylinders in blue phases, the bending 
energy of the columns (energy ruled by the Prank constant ^3) becomes prohibitive for large 
radii. Which distance I from the axis is difficult to guess at this stage; in analogy with BP's, 
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one could take I = zu, which yields r^ax i^)"^ — vj\/w^^^^}? (^max(A) is not the distance of a 
helix to the axis, but the span of the helicoid; see equation 5 and appendix). On the other 
hand the observations reported in [12] indicate that the domain of existence of the HR is not 
limited by a cylinder; taking then rmax — ^ (independent of A), one gets Z(A) = ^Jxu'^ -\- 
Equation 12 is then replaced by 

< r < Tmax = Ct7; arctan — <Lp < arctan — — . (loj 

A A 

(b)- A solution with positive Gaussian curvature has an obvious disadvantage: because 
the columns have a finite width, a whole number of columns does not fill the interior of a 
layer without void or overlapping; and in case it succeeds to fill the interior of some layer 
without void or overlapping, it does not in the next one. Therefore the columnar packing 
necessarily exhibits edge dislocations, whose energy will be estimated later on to be much 
larger than the energy of the helical geometry. Contrariwise, the columns have space to 
move locally along a lateral direction if the layers have Gaussian negative curvature, and in 
fact to take in the layers a distance of equilibrium that can be different from the distance 
from a layer to the next, without necessitating edge dislocations. Thus, parallel helicoids 
are preferred to nested cylinders, although the distribution of the columns is rotationally 
invariant (in a continuous analytical description). 



III. THE FREE ENERGY 

As already stated, the ribbons, according to [12], are characteristic of metastable phases 
B7*; which in addition to Sm and D orders, display helicity (H order), either left or right. 
Thus the finite pitch p {— 2Tr/q) of the columns is a material constant (which might depend 
on the temperature). Assuming this to be true, we show now, in the framework of a very 
simple model, that wq = — 1, i.e. b = —p. 

The free energy of a HR contains a part relative to the columns, with density Cch, total 
energy Ech, which can be written as the free energy of a cholesteric phase oriented along it, 
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and a part relative to the layers, layer energy density eiay, total energy Eiay, which can be 
written as the free energy of a smectic phase. We shall assume that the phenomenological 
free energy density is a sum of these two contributions, without any term of interaction. All 
the interaction is indeed contained in the isometric geometry requirements: the columns are 
equidistant and the layers are parallel. For the same reason, we do not introduce any energy 
of elastic (positional) deformation. 

The helicity of a HR is described by two pseudoscalars 

(a) - the Burgers vector b of the dislocation, b — ndo, n G Z, do is the smectic repeat 
distance. We take b positive if the dislocation is right screw, negative in opposite case; 
zu — 6/27r. 

(b) - the pitch of the the twisted columns. This pitch is perpendicular to the layers. We 
take the pitch positive if the helicity of the columns is right handed, negative in opposite 
case; q — 27r/p. 

The 'cholesteric part' of the energy reads 



1 1 
^ch — 7.^2 • curllt + g) + -K^ {'rt x curllty 



2K2 - 



qK2 + -^Ks 

4:ZU 



(16) 



Observe tat if ip'^ > 0, one expects q > 0, because the twist q contributes to a reduction 
of the energy density. Therefore, according to our discussion of the signs above, the inverse 
pitch g, and the Burgers vector b (which is the pitch of the helicoids) are of opposite signs. 
At this level of calculation we do not specify the origin of q when assigning a cholesteric 
energy density to the director field of the B7* phase. In other words, here q is not necessarily 
a material constant (although such a possibility is not excluded) but rather an empirical 
model parameter that reflects the (metastable) helicity of the columns in the B7* HRs. The 
result given below shows that q is intimately related to the Burgers vector b. 

The total 'cholesteric' energy per unit length along the axis is 

Ech = -^l 1 1 echdXdS, dS = -^NdrdO; (17) 
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dS is the element of area of H( A) . When integrating along one pitch of the hehcoids, 9 
runs over the interval [0, 27r[; it contributes by a factor 2% to the integration, since no term 
depends on in 9 the integrand. The rest of the region of integration is a variable of the 
problem, which is adjusted by the competition between layer curvature effects (discussed 
below) and cholesteric contributions. In this section, we just assume that it scales like 
zu"^. This region is necessarily bound by the two sheets of the focal surface F (the loci of 
the centers of curvature of the helicoids) on which tan2(/? = 1. Starting from the central 
hehcoid H(A = 0), the first hehcoids for which this equahty is obtained are H(A = w) and 
H(A = —zu), at r = (i.e. on the straight generators of H(A = 0), at a distance ±zu of the 
central axis). For the time being, and in order to avoid a difficult analytical calculation, we 
write ~ (ech) with {(p'^) = -|ti7, {ip'^) = \zu^, i.e. 

11 r 1 " 

E^h^{ech)^^^^K2 + -Ks + K2qzu 1 + -qw (18) 

If this expression were the only energy term, one would get, by minimizing with respect 
to 77 = zuq: 

wq = -1 (19) 

The sign of the vo and of the q are opposite. It is easy to check that {ech) is a minimum 
for this value of va {i.e. ^-^^ > 0). The contribution of the layer curvature energy does 
not change much these results, as we show now. 

The layer curvature of the energy density reads: 

eiay = 2k {Hf + K„G (20) 

K and K}) are referring to one layer. The total energy can be written: 

f f f dSdX , , 

^lay = J J J ^lay^, (21) 

do being the layer thickness. We discuss the k and Kb term separately. 
Mean curvature. The set of stacked layers is built upon a central layer that is a mini- 
mal surface. This situation is most probably generic, because it yields such a small mean 
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curvature energy. We shall therefore neglect the 'regular' part of this energy, keeping only 
the 'singular' part, that originates in the presence of the k = 1/2 singularity at r = 0, 
X = zu. The disclination length is b\/2 per turn. We thus expect an energy of the order of: 
^^^'^ In ^ per unit length of helical axis. 

Gaussian curvature. Assume that in each layer the boundaries of the domain of existence 
of the phase are along two columns ri, r2. According to the Gauss-Bonnet theorem we have 
/ / GdS — n — J GgdS, with Ug the geodesic curvature. The geodesic curvature vanishes 
for a geodesic, and takes the value = ±^ cosip = ±^ simp on a helix (a column), see 
appendix. The length of a column parameterized by r is 1/ sin a per unit length of helical 
axis, a (r) being the angle of the tangent it to the column with the plane orthogonal to 
the HR axis; sin a — riz = ^ sin (p, which is positive according to our conventions. We thus 
expect an energy of the order of Kb f ~ + per unit length of helical axis for the 
layer H(A). Because N'^ — r'^ + w"^ does not depend on A, the contribution is the same for 
each layer, and the expression for one layer has to be multiplied by the number of layers 
w jd^. Take ri = 0, r2 = cc^, one gets | — + = Kb'K^^^^. Observe however that 
the contribution of the disclination is certainly much larger than the contribution of the 
Gaussian curvature, the ratio being of the order of Kvo/K^dQ, where ^ =1, and ^ >> 1. 
We therefore neglect the Kb contribution and eventually get: 

£„^^,„-=I^,„- (22) 

do do 4do do 

This exspression is of the same order of magnitude as, althoug presumably smaller than 
(see below, after Eq.25) Ech, the ratio being of the order of Ech/Eiay = K2do/K. 

We now minimize the total energy E = Eiay + E^h with respect to the variable rj = qzu. 
One gets: 

dE ttV2k 1 , d'^E nV^K 1 

-^ = ^ + ir2 1 + 77 ; -^ = -^ J + 23 

OTj 4do Tj OTj^ 4do T]^ 

the analysis of the stability of the solution (which obeys ^ — 0), starts with the condition 
> 0, in which inequahty one inserts the value of r] for which §^ = 0. After some algebra, 
one finds one stable solution only 



that gives back equation for k = 0. Notice that the existence of such a solution requires 
that 

1 - TT^-^ > (25) 

i.e. that n/dg be smaller than K2. This seems to be generically the case. One can indeed 
argue as follows. Due to the anisotropy of the splay modulus, the layers most certainly fold 
more easily about the columns than perpendicularly to the columns. Call ap the curvature of 
a planar section of H(A) perpendicular to the columns {ap measures the folding of the layer 
about the columns), Uc the curvature of a section parallel to the columns {ap + ac — 2iJ), 
one finds 

a„ 1 + n 

^ = -1 + 2--^, t'>l, 



SO that 



> 3. Therefore the layers are essentially folded about the columns, and the splay 
modulus at work is small. Furthermore K2, which is the twist constant of long (in principle 
infinite) columns, is expected to be large. Therefore the inequality 25 is certainly satisfied, 
and wq of Eq.24 is a stable solution. 

It remains to show that the above solution is much less costly than a solution with circular 
cylindrical layers, which would necessarily contain edge dislocations of the columns (with 
Burgers vectors dc, which is the columnar repeat distance. The density of such dislocations 
would be pe ~ 1/Rdc, at a distance R from the cylinders axis. Therefore the number 
of dislocations in a cylinder of radius zu (which scales like their length per unit length of 
cyUnder) would be pe ~ / RdRd9/rdc — ^i^w/dc- The associated energy is of the order 
of per unit length of dislocation, i.e., E^isi ~ 2'KBdc'W per unit length of cyUnder. 
Comparing with the energy of the HR calculated above, one gets Ehr/ E^isi ~ In 
this expression, K2/B is a squared microscopic length, of the same order as d^. Therefore 
Ehr/ Edisi is the ratio of a microscopic length over a macroscopic length, and is very small. 
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IV. DISCUSSION 



Isometric textures are the rule in liquid crystalline phases as soon as there is one or two 
dimensions of periodicity [4]: focal conic domains when the order parameter is layer- like 
({Sm}, one dimension of periodicity), developable domains when it is columnar ({D}, two 
dimensions of periodicity). In both cases the associated distorted textures are metastable 
and carry singularities. The situation is more involved and really fascinating when is added 
to the periodicities a trend to helicity, {H}. Twist grain boundary phases (TGB) [24] and 
Blue Smectic phases (BPSm) [25] offer remarkable examples of the {Sm+H} mixed case; in 
both cases the hehcity axis is in the layer and affects the molecules. Then, one gets frustrated 
phases, regular, periodic, with thermodynamically stable singularities (dislocations, mostly). 
{H+D} frustrated phases have been suggested in [26]; they could display domains of parallel 
columns separated by twist walls, possibly observed in [27]. 

The {Sm+D} order parameter is such that the layers are modulated in thickness and fit 
closely along the modulations. We have demonstrated that helical ribbons are manifestations 
of the isometric principle for this order parameter. Those HRs have a finite width and do not 
transform in another texture when they come into close contact. For instance, they do not 
extend laterally into quasi-planar layers, after the manner of the giant screw dislocations 
in SmAs, alluded to above. This extension would require a slight breaking of isometry, 
and the columns would no longer be twisted, just bent. It is probably this latter character 
that forbids such an extension. In other words the helicity seems to belong to the order 
parameter, which is therefore {Sm+D+H}. But there is no way to construct in Euchdean 
space an infinite crystal with helicity showing the same close packing as in the (Sm+D) 
case; the HRs can then be said to be frustrated textures. The {Sm+D+H} phase is most 
probably crystal-definable in a curved space, but we do not have investigated this question. 
Note however that the isometric principle could be directly applied to this curved crystal, 
and would then yield the HRs. 

It is worth mentioning another character which differentiates SmAs giant screw dislo- 



20 



cations and HRs. While the HR geometry defined in section II is experimentally strictly 
obeyed (we have already mentioned the rigidity of these textures) , the inner geometry (inside 
the helical disclinations) of a Sm screw dislocation is very flexible and susceptible of many 
deformations. In fact, the giant screw dislocation is only an approached isometric texture 
for a Sm phase (the true isometric textures are focal conic domains) , and can therefore suffer 
large fluctuations, in particular due to the presence of edge and screw dislocations of small 
Burgers vectors ( [28], chapter 5). 

We have mentioned the analogy between a HR and a double-twist BP cylinder. In both 
cases the columns (in the BP case what we call columns are the envelops of the director 
field) are double-twisted and rotationally invariant about the axis direction. Because BP is 
a true 3D liquid, there is no obstacle to merging when two double-twist cylinders come into 
contact; the director field deforms without much energy expense between cylinders, and the 
cyhnder geometry fills space [29]. Now, as already stated, a layer hehcal texture is preferable 
to a layer cylindrical texture for a HR. There are no layers in a BP, and it is the custom to 
desribe this double-twist texture as if the columns are partitioned into cylinders. It is true 
that isometry requirements are weak in BP, as long as the molecules forming the columns 
are short compared to the pitch. However, it is believed that BPs are favored with long 
molecules [30,23], and double- twist cylinders are also present in certain DNA solutions [31]. 
One might therefore wonder whether a partition into parallel helicoids (compared to nested 
cylinders, this requires a small displacement of the columns) might not be closer to reality 
in some BPs. 

A final, brief, remark relating to the HR width limitation. Let us notice that the region 
of negative Gaussian curvature is bound by the focal surface of the stacking, whose two 
sheets obey the equation t = ±1, see section II, C. This is a first reason why the width 
stops short of the focal surface. The helical disclinations belong to the focal surfaces, which 
extend outwards. A second reason relates to the twist of the columns: notice, equation 11, 
that it -curllt ~ 1/A for A large, i.e. the twist vanishes for the helicoids that are far enough 
from the axis; consequently the twist energy becomes large for A large, if, as assumed, p 
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is a material constant (metastability of the B7* phase); at the same time the bend energy 
decreases. Observe that these behaviors of the twist and the bend energy for A large are 
both opposite to their behaviors for the cylinders of double twist. But in both cases one 
expects that the competition between these two energy terms fixes the width of the texture 
at some finite value which scales with the pitch. 
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V. APPENDIX 

a) width of a HR. 

Let tan'^ = Xvo/Nr. Eq.5 can be written 

H{\) = {x^lcos{9-i;), y^lsin{e-il;),z^w{e-'ijj)+r](l)}, (26) 

where P — r'^ + is the distance to the z-axis; rj (l) — wt/j + ^ tan'^. 

Assume now that rmax = ^, one gets: 



2 , , / „ A2 zu /, N , A 

tan ^ = — — ; /max = \h^^ + tt = V (w) = zu^ + 

2zu^ V 2 cosV' V2 

Therefore the equation of the boundary of the nucleus in a plane z ^constant (here this 
constant is taken equal to zero) is, in polar coordinates: 

h) geodesic curvature. 

Let M be a 2D manifold, M be a point on M moving along a curve C, s the abscissa 
along C, u the normal in M to M, we have, according to well-known formulae of surface 
theory: 

ds^ ^ \ ds j ds^ 
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The normal curvature Gi, is the curvatureof the section normal to M at M; its sign 
depends on the sign of . The geodesic curvature (jg measures the twist of the curve in 
the plane tangent to M in M. It vanishes for C a geodesic line. Its sign depends on the 
orientation of the curve. 

These quantities are easily calculated for a helix on the hehcoid H(A) . We have — n . 
Thus, follows from the expression of rf by differentiating at constant r. The calculation 
yields: 



Hence: 



^ sm^ {-%, n^, 0} (= curllt x it) 



1 r 

(Jv — —\ COS^ </? — \ — COS (/?. 

A Xzu 
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Figure Captions 

Figure l.Edge dislocations in a bent crystal: (a) at random, (b) polygonized; (c) curved 
layers in a SmA, infinitesimally small edge dislocations, isometric deformation. 

Figure 2. The banana-like chemical configuration of D14F3. 

Figure 3. A slim hehcal ribbon (S-HR) observed under the optical microscope. 

Figure 4. Experimental data relating to the (slim) Helical Ribbons: (a)- distribution of 
the widths; note that the experimental accuracy is ~ lyum, so that each box is 2/im wide; 
(b)- ratio h/l for the same experimental data; this ratio does not vary much when the width 
deviates from its more probable value. 

Figure 5. Geometry of a helical ribbon: (a)- nested helicoids stacked parallel to a central 
ruled helicoid; the domain width is hmited |r| < &/27r (see Eq. 5); the 'holes' in the 
figure result from the small number of helicoids chosen parallel to the central HR; (b)- 
the A = (central helicoids) and A = 0.5 helicoids, with the modulations (the columns) 
along hehcal lines; |r| < 6/27r. In each figure the two helices are the singularities of the 
stacking (dischnations k — 1/2). The Burgers vector is assumed positive; the hehcoids are 
right-handed. 
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